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Background srw 

Last Year's AHS Forum (V22 TRAM Rotor)* 

•Demonstrated improved prediction of figure of merit (FM) 

❖ Improved blade grid resolution with high-order spatial differences 

❖ Spalart-Allmaras turbulence mode with the DES length scale 
•Demonstrated Cartesian off-body (OB) adaptive mesh refinement (AMR) 

❖ No affect on FM and blade loads 

❖ Resolved vortex wake rich in turbulent flow physics 

❖ Vortex worms: Entrainment of shear layers by vortices 



*Chaderjian, N. M. and Buning, P. G., "High Resolution Navier-Stokes Simulation 
of Rotor Wakes, AHS 67 th Annual Forum, Virginia Beach, VA, May 3-5, 2011. 
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Objectives 
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•Demonstrate on the UH-60 flexible rotor in forward flight and hover 

❖ OVERFLOW'S high-order spatial differencing, AMR process, and SA- 
DES turbulence model 

❖ Convergence of the dual-time algorithm 

❖ Convergence of loose coupling process with CAMRAD-II and AMR 

❖ Explore details of the rotor wake in forward flight and hover 
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Numerical Approach For This Study SRW 

OVERFLOW Version 2.2 

•Time-dependent Navier-Stokes equations solved everywhere (body & off-body) 

•Overset body-fitted curvilinear grids rotate through fixed off-body Cartesian grids 

•2 nd -order time accuracy using dual time-stepping (pseudo-time with sub- 
iterations) 

❖ Approximately factored Pulliam-Chaussee diagonal algorithm 
•5 th -order convective differences, 2 nd -order grid metrics/viscous terms 
•Dynamic adaptive mesh refinement (AMR) for the off-body (OB) Cartesian grids 
•OVERFLOW is loosely coupled with CAMRAD-II every % revolution (4-bladed rotor) 

❖ CFD delta-loads replace CAMRAD-ll's simplified aerodynamics model 

❖ CAMRAD-II provides OVERFLOW with blade deflections and trim conditions 
•Many other options and features 
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Spalart-Allmaras Turbulence Model 
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"d " is the distance to the nearest wall 

• It should represent the largest turbulent eddy size 

• It's too large in the wake, so TEV is too large 

* Production, but no dissipation (out of balance) 

• In the wake it has simply become a geometric 
parameter 
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Spalart-Allmaras Detached Eddy Simulation SRW 
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DES provides a more realistic turbulent length scale 

d = min(d,C DES A) 

A = max (Ax, Ay, Az) 

• Rational way to reduce the TEV, even when using RANS grid spacing 

• Viewed as a hybrid model 

❖ RANS in boundary layer 

❖ LES outside of boundary layer 

• Turbulent structures are resolved if the mesh is refined 

❖ Unique application of DES to rotorcraft wakes 

❖ Vortices are not diffused 

❖ Blade flow separation still determined by RANS model 

• Use delayed detached eddy simulation (DDES) as a precaution 
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Overset Grid System - UH-60 Rotor 
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Modified 

Shank 


Component 

Dimensions 

Grid Points 

Rotor Blade 

241x291x63 

4,418,253 

Inboard Tip 

137x69x63 

595,539 

Outboard Tip 

122x49x63 

376,614 


One Blade Total 

5,390,406 

Hub Body 

93x301x54 

1,511,622 

Hub Top 

21x21x54 

23,814 

Hub Top 

21x21x54 

23,814 


Hub Total 

1,559,250 
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High-Speed Level Flight 

Low BVI 
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I — ► r/R Flight-Test Measurements 




Flight Parameters: C8534 
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Establishing Dual Time-Step Asymptotic Range 
2 nd -Order Accuracy: C8534 
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Normal Force Waveform 



•RMS difference formed from highly converged result 
•Asymptotic range begins at Nj=35, At=%° 

❖ Typical choice: Nj=15, At=}T 
•Differences are quantitatively larger than they appear 
•Chose Nj=40, At=%° 

❖ 2.6-order sub-iteration drop is larger than the 
2-order rule-of-thumb 


Convergence with Timestep 






SRW 



AMR Grid Systems: C8534 

AMR efficiently captures wake without grid refinement 
❖ 80 million grid points with uniform 10% c tp grid 
AMR allows for efficient/automatic grid refinement 



AMRO: A=10% c tiD , 960 grids, 61 million grid points 
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OVERFLOW/CAMRAD-II Loose Coupling Convergence 

Control Angle Convergence 
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•AMRO: Run for 4.5 revs 

•AMR2: Run for 5.0 revs (Started from AMRO, 3 revs) 

•Both AMRO and AMR2 are converged at coupling iteration 12 
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RMS Difference (A Iteration) 



OVERFLOW/CAMRAD-II Loose Coupling Convergence 

RMS, Normal Force and Pitching Moment Coefficients 


SRW 


•RMS difference between successive iterations (ANJ 
•Both AMRO and AMR2 RMS values are converged at cycle 12 

•Both AMRO and AMR2 loads show little change past 12 th cycle 

Convergence of Normal Force 
and Pitching Moment Coefficients 


AMRO AMR2 






OVERFLOW/CAM RAD-II Pitching Moment Coefficient 

Flight Counter C8534, Mean Removed 

AMRO and AMR2 results virtually identical 
❖ This result is expected due to the lack of BVI 
Good agreement of RMS differences with flight-test data* (to 2.5%) 
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*Gloria K. Yamauchi and Larry A. Young, "A Status of NASA 
Rotorcraft Research," NASA/TP-2009-215369, September, 2009. 
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OVERFLOW/CAM RAD-II Normal Force Coefficient 

Flight Counter C8534, Mean Removed 
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• AMRO and AMR2 results virtually identical 

❖ This result is expected due to the lack of BVI 

• Good agreement of RMS differences with flight-test data (to 2.1%) 



*Gloria K. Yamauchi and Larry A. Young, "A Status of NASA 
Rotorcraft Research," NASA/TP-2009-215369, September, 2009. 




Flow Visualization Using the Q-criterion 

AMR2: A=10%, 5% and 2.5% 
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Flight Counter C8534 
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Flow Visualization Showing AMR Tracking Vortices 

Rendered on Horizontal and Vertical Planes 
AMR2: A=10%, 5% and 2.5% 
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Time-Dependent Surface Flow 

Texture Mapping Colored by Pressure 
AMR2: A=10%, 5% and 2.5% 
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OVERFLOW/CAMRAD-II UH-60 Rotor in Hover 
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AMRO: A=10% 

1,800 grids, 78 million grid points 



AMR1: A=10% and 5% 

7,700 grids, 302 million grid points 



17 




SRW 



OVERFLOW/CAMRAD-II UH-60 Rotor in Hover 

AMRO and AMR1 FM in close agreement with each other 

❖ 2% difference from UTRC test 

❖ Experimental uncertainty unknown 

AMRO and AMR1 sectional thrust coefficient in good 
agreement with UTRC test 
CFD accuracy has improved 

❖ Further validation requires test data uncertainties 
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Flow Visualization Showing AMR Tracking Vortices 

AMRl: A=10% and 5%^ 

7,700 grids, 302 million grid points 
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Cut-Away View of Vortices and Worms 

AMRl: A=10% and 5% 

7,700 grids, 302 million grid points 
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Pressure Iso-surface for the UH-60 in Hover 

AMRl: A=10% and 5% 

7,700 grids, 302 million grid points 

•Pressure iso-surface similar to viewing natural condensation of tip vortices 

❖ You can not see any worms (Typically 0.5%-l% of tip vortex strength) 

❖ Worms are evident by bands of vorticity on the iso-surfaces 
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Computer Run Times 
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•OVERFLOW was executed on NASA's Pleiades supercomputer system using 
Intel 2.93GHz Westmere nodes and MPI 

•Each Westmere node has 12 cores and 24GB of memory 

•Wall clock timings are not optimal, but reflect what is possible given 
computer resources and deadlines 


Case 

Grid Points 
(millions) 

Number of 
Cores 

Wall Clock 
Time (hrs/rev) 

Relative 

Processing 

Rate 

C8534 

AMRO 

61 

1,536 

5.4 

0.13 

C8534 

AMR2 754 3 ' 072 238 ° J1 

Hover 

AMRO 

78 

1,536 

5.8 

0.15 


Hover 

AMR1 


302 


4,608 


10.1 


1.00 
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Conclusions 
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• Time-dependent Navier-Stokes simulations have been carried out by loosely coupling 
OVERFLOW with CAMRAD-II. 

• Time accuracy of the Navier-Stokes equations was established within the asymptotic 
range, but required more sub-iterations/convergence than normally employed. 

• OVERFLOW/CAMRAD-II loose coupling rapidly converged, even with AMR. 

• AMR can be used to find vortex wakes without an apriori knowledge of the wake 
shape or location, allowing solutions with fewer grid points. 

❖ Useful for baseline or highly resolved wake simulation. 

• Computed RMS differences of M 2 C N and M 2 C m for forward flight agreed with flight 
test data* to 2.1% and 2.5% respectively. 

• Computed FM agreed with experiment to 2% 

❖ More experimental uncertainty data is needed for further CFD validation. 

• Complex turbulent wakes interacting with vortices formed vortical worms, similar to 
last year's V- 22 simulations. 

• AMR had little affect on the rotor loads for forward flight due to weak BVI. 

*Gloria K. Yamauchi and Larry A. Young, "A Status of NASA 
Rotorcraft Research," NASA/TP-2009-215369, September, 2009. 
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UH-60 rotor/hub in hover 
AMR1: A=10% and 5% 

7,700 grids, 302 million grid points 
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